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ABSTRACT 


The  dynamics  of  phase  locked  loops  have  been  studied  by  many 
investigators  for  the  cases  of  either  a  pure  sinusoidal  signal  or 
for  a  pseudo  noise  (PN)  sequence.  This  dissertation  presents  an 
analysis  of  a  phase  locked  system  consisting  of  a  sinusoidal  carrier 
modulated  with  a  PN  sequence.  The  effect  of  the  modulation  index  of 
the  phase  shift  modulation  and  the  design  of  a  loop  filter  for  coherent 
detection  of  a  PN  sequence  are  analysed.  Additionally,  the  error 
probability  for  each  digit,  \<ihen   a  likelihood  detector  is  used  in 
front  of  the  delay  lock  discriminator  for  a  PN  sequence,  is  also 
converted  to  the  delay  error  of  the  phase  lock  loop. 
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CHAPTER  I 
INTRODUCTION 

A  correlation  type  detector  or  a  matched  filter  is  known  to  be 
optimum  for  detecting  a  signal  in  additive  noise  under  a  vide  variety 
of  conditions.  The  noise  power  is  proportional  to  the  effective  band- 
width for  the  case  of  white  Gaussian  noise.  We  can  improve  the  s/N  ratio 
by  narrowing  the  effective  bandwidth  if  it  does  not  restrict  the  infor- 
mation contained  in  the  signal. 


s(t)s(t+T) 
s(t)+n(t)    , ,  -|-n(t)s(t  +  T) 


Multiplier 


Low  pass 
Filter 


y(t) 


s(t+T)| 
Fig.  1-1   A  Correlation  Type  Detector 

Consider  the  detector  in  Fig.  1-1.  s(t)  is  the  incoming  signal  and 
n(t)  is  additive  white  Gaussian  noise.  If  we  know  the  signal  except  for 
its  phase,  that  is,  the  phase  of  the  signal  is  the  only  information 
contained  in  the  signal,  we  can  generate  sCt+T)  to  multiply  with  the 
incoming  signal.  The  low  pass  filter  works  like  an  averager  in  some  sense 
as  described  later  in  section  4-3,  and 

y(t,.T)=YJs(t)s(t  +  r)dt+:j:jn(t)s(t-f-T)dt. 
If  s(t)  and  n(t)  are  assumed  to  be  stationary,  y(t,,T)  becomes  a  function 
of  only T for  long  T,  and  y(t)  becomes  approximately  the  sum  of  the  auto- 
correlation function  of  s(t)  "R^(t)"  and  the  crosscorrelation  function 
of  s(t)  and  n(t)  "R^p(T)". 

y(r)^  R^(T)+R^„(r) 

If  we  can  find  a  signal  such  that  its  autocorrelation  function  has 
a  distinguished  difference  between  the  values  of  T=0  and  ff  0,  we  will 
be  able  to  utilize  it  for  phase  detection, 
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A  sine  function  {z   A  cos  a)c.t  has  the  autocorrelation  function  of 
R^iX)  — A^oosuJcXt   and  the  phase  locked  loop  for  a  sinusoidal  function, 
like  that  shown  in  Fig,  1-2,  has  been  studied  by  many  investigators. 
Among  them,  R.  Jaffe  and  E,  Rechtin  showed  an  optimum  filter  design  for 
a  linearly  approximated  loop  [5  J  t  A,  J.  Viterbi  showed  the  transient 
response  for  Ui   having  a  second  order  polynomial  in  the  nonlinear 
system, [lO]  W,  F,   Roland  defined  the  threshold  criterion  by  the  minimum 
of  the  cycles  of  slip  and  shovjed  its  utility  through  experimental  results 
in  his  thesis  for  the  degree  of  M.  S.  [?] 

V2Acos(u;ct  ■^a/)+n(t)     ACk„sin(9r0^)^"\/2Ck„n(t)sin(M;ct  t  ^J 


v/2"Csin(M;ct+92) 


Low  Pass 
Filter 


Voltage  Controlled 
Oscillator 


Fig.  1-2  Phase  Locked  Loop  for  Sinusoid 
The  binary  signal  has  only  two  possibilities  for  each  digit,  ONE  or 
ZERO,  or  -1  or  -f-1,  but  sequences  of  binary  digits  show  many  interesting 
and  ccanplex  characteristics.  Some  sequences  represent  codes  for  which 
maximum  likelihood  detection  is  possible.  They  are  classified  by  the 
crosscorrelation  function  in  each  code  set;  orthogonal,  biorthogonal  and 
simplex*  [2 J 

There  are  also  some  special  group  of  binary  sequences  especially 
adapted  for  phase  detection;  the  PN  sequence,  especially  the  maximum- 
length  linear  shift  register  sequence  (m-sequence),  will  be  explained  a 
little  more  precisely  in  Chapter  2.  R.  H.  Barker  defined  "Group  Synchro- 
nization" so  as  to  pin-point  effectively  the  origin  in  time  for  a  binary 
sequence,  and  showed  a  matched  filter,  and  also  the  characteristics  and 
some  patterns  of  this  sequence,  [ij 

The  complementary  binary  sequences  studied  by  S.  Jauregui,  Jr.  used 
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the  same  kind  of  matched  filters  and  when  used  in  a  differential  pair 
these  eliminate  small  ambiguities  in  the  output.  ^6]  Both  sequences  can 
be  used  in  the  form  of  discrete  pulse  sequences. 

PN  sequences  are  adequate  for  the  use  of  both  continuously  a.nd 
periodically  incoming  signals.  J.  J.  Spilker  proposed  a  delay  lock  dis- 
criminator shown  in  Fig.  1-3 »  using  the  combination  of  two  autocorrelation 
functions  having  two  digit  interval  difference  T.  Ls]  J.  C,  Springett 
showed  a  combination  of  clock  and  PN  sequence  to  produce  a  controlling 
voltage  proportional  to  the  phase  difference.  [9]  M.  F.  Easterling 
introduced  a  tracking  system  using  a  combination  of  clock  and  several  codes 
and  produced  synchonization,  code  by  code,  separately.  [3]  In  this  system, 
As(tfT)-(-n(t) 


IxDop  Filter 


|C/)(t-^+t) 


Voltage  Controlled 
Clock 


Linear  Feedback  Shift 
Register 

Fig.  1-3   Delay  Lock  Discriminator  for  Binary  Signal 

only  the  clock  has  a  phase  locked  loop.  The  synchronization  of  the  other 

codes  is  evaluated  by  a  correlation  indicator.  For  the  system  which  needs 

long  distance  tracking  and  does  not  require  quick  reaction  for  the  variation 

of  the  targets  dynamics,  this  will  be  all  right,  but  if  we  need  autanatic 

synchronization  and  very  rapid  acquisition,  this  will  not  be  adequate. 

Since  the  second  one  has  many  weak  ambiguous  points  between  distin- 
guished phase  lock  points,  v/e  prefer  Spilker 's  delay  lock  discriminator 
rather  than  Springett 's. 

The  periodic  sequence  is  preferable  for  phase  lock  loops  to  the  non- 
periodic  because  of  the  possibility  of  continuous  control. 

This  paper  treats  the  phase  lock  system  for  a  carrier  phase  shift 
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modulated  by  a  PN  sequence,  as  shown  in  Fig.  l-4o  s(t)  is  a  maximum- 
length  linear  shift  register  sequence  with  the  number  of  digits  in  a 
period  M=2  -1,  and  the  duration  of  a  digit  d.  n(t)  is  white  Gaussian 
noise  and  has  the  following  narrow  band  approximation 
n(t)=  n^(t)cos«(t  +  no(t)siniJ^ 


Figo  1-^  A  Phase  Lock  System  for  a  Signal  Modulated 
by  a  PN  Sequence 

Many  papers  have  analy^sd  each  separate  stage  presented  so  far,  but 
not  the  combined  system.  There  are  two  problems  to  be  solved  for  this 
system;  what  is  the  effect  of  9  on  the  phase  lock  performance  of  this  system, 
and  if  FiC'^)  and  F^(u>)   have  the  same  type  of  low  pass  filter  as  that  derived 
by  Jaffe  and  Rechtin,,  hoif  can  we  deterraine  adequate  bandwidths  for  these 
filters?  9  determines  the  ratio  of  the  I  component  and  Q  component  of 
incoming  signals,  and  also  the  necessary  bandwidth  for  the  input  of  this 
system.  If  0  is  large,  the  output  modulating  signal  will  be  large,  but  due 
to  the  large  phase  error  of  the  demodulating  stage,  the  effective  amplitude 
will  decrease.  The  loop  bandwidth  of  each  stage  is  related  to  the  tran- 
sient response  for  tracking  and  the  required  reduction  in  the  noise  output. 

This  report  will  summarize  those  relations  determined  by  other  authors 
for  each  stage,  and  will  further  derive  the  threshold  s/n  ratio  and  the 
performance  of  this  system.  Chapter  2  considers  the  characteristics  of 
m-sequences.  Chapter  3  shows  the  demodulation  of  an  incoming  modulated 
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carrier  and  the  output  of  this  stage. 

The  delay  lock  discriminator  for  an  m- sequence  will  be  analyzed  in 
Chapter  4  with  the  effect  of  a  llmiter  in  front  of  this  stage  included 
in  the  performance. 

Chapter  5  summarises  the  system  results  and  illustrates  by  an  example 
the  conclusion  of  the  next  chapter.  This  illustration  will  not  show  any 
particularly  interesting  value  former  other  parameters.  It  is  simply  an 
example  of  the  evaluation  of  system  parameters  and  of  design  decision, 
but  it  is  hoped  that  this  will  give  a  rough  idea  of  the  overall  performance 
of  such  systems  and  will  give  some  understanding  of  the  role  played  l^  each 
parameter  in  the  design  of  a  system  employing  PN  modulation  of  a  sinusoidal 
carrier. 
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CHAPTER  II 
MAXB'iTOI-LENGTH  LINEAR  SHIFT  REGISTER  SEQUENCE 

Tlio  autocorrelation  function  of  a  binary  sequence  of  length  M, 
s(i)  =  (a,  ,a^,  0..  sa^j,  v/here  a^  =  ±l»  is  defined  as 
&(3)=pit^ia.i,j         itO  5  modulo  M 


M 

I 

(=1 

If   £.(0)=1 


^(j)=a  +  l    (l<j<M) 
then  s(i)  is  called  a  two-level  sequence,  and  if,  moreover,  ^(j)=--ir 
for  j^O  (modulo  M) ,  the  sequence  is  sometimes  called  a  pseudo-noise 
(PN)  sequence.  Exactly,  a  PN  sequence  has  to  have  the  following  three 
"randomness  properties".  [4J 

R=l  (the  balance  property)  In  each  period  of  the  sequence  the  number  of 
ONE'S  differs  from  the  number  of  ZERO'S  by  at  most  1. 
R-2  (the  run  property)  Among  the  inins  of  ONE's  and  of  ZERO'S  in  each 
period,  one-half  the  inins  of  each  kind  are  of  length  one,  one  fourth  of 
each  kind  are  of  length  two,  one  eighth  are  of  length  three,  and  so  on 
as  these  fractions  give  meaningful  numbers  of  runs. 

R-3  (the  correlation  property)  If  a  period  of  the  sequence  is  compared, 
term  1:^  term,  with  any  cyclic  shift  of  itself,  the  number  of  agreements 
differs  from  the  number  of  disagreements  by  at  most  1, 

There  are  four  types  of  PN  sequences,  corresponding  to  the  length 
M»  [2] 

1)  M=2"  -1  (m=sequences) 

2)  M=^4t-1  is  prime  (quadratic  residue  or  Legendre 

sequence) 

3)  M=4t-1  =  4x^+27  is  prime       (Hall  sequence) 

4)  M=p(p+2)  where  both  p.,  p+2  are  prime  (twin  prime  sequence) 
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Among  them  the  maximum-length  linear  shift  register  sequence  (m-sequence) 
has  the  delay-and-add  property  and  is  generated  by  a  linear  shift  register. 
This  is  also  convenient  to  generate  in  the  phase  lock  loop,  therefore  in 
this  report  the  m-sequence  is  used  as  a  PN  sequence, 
clock 


r—j 


Xj 

r 


X3 


1 


1  ^X  I  fX^f  Xj  t  »  »  »  ,Xn^ 


f 


->  Out 


g  ^X I  ,X  2 1  Xj  » •  •  •  »x^ ) 


clock 


il       kl       il 


■^Out 


f(x,  ,X2,  x,,"',:;^ 


Fig,  2-1  n  Stage  Linear  Shift  Register 
In  the  n  stage  linear  shift  register  shown  in  Fig.  2-1,  the  output 
sequence  is  always  periodic  with  a  period  p  ^  M  =  2  -1 ,  and  when  period  is 
M,  that  sequence  is  called  a  maximum-length  linear  shift  register  sequence. 
This  m-sequence  has  the  delay-and-add  property,  that  is,  if  the  sequence 
is  added  term-by-term  modulo  2  to  a  cyclic  shift  of  itself,  the  resulting 
sequence  is  another  cyclic  shift  of  the  original  sequence.  This  property 
is  possessed  only  by  the  binary  sequence  {a,  ,ai,.. .  ,a„}  (ai=l  or  0)  which 
satisfies  a  linear  recursion  relation 

for  all  k 
k-i  :  (modulo  M) 
Ci=  1  or  0 

summation  is  modulo  2 
and  M  =  2*^  -1.  (.2]  m-sequence  satisfies  this. 

Consider  an  m-sequence  s(i)s={a,,  a^,...,a„|  .  The  autocorrelation 
function  of  this  sequence  ^(k)  is 

A(k)=^^ai®a.+j^  0,  sum  of  modulo  2 


Using  the  delay-and-add  property, 

^^M  =  -p\Z^i*]   =h{(#  o^  0NE's)x(»1)+(#  of  ZERO's)x(l)] 


I*/ 


=  it^x(.i).^x(.i)/ 

—  =  I  for  kf  0  (modulo  M) 

^x.(0)=j^iaiea,  =  ji  1.1=1 
(If  we  replace  ONE  and  ZERO  with  =1  and  +1,  we  can  use  multiple  instead 
of  0 )o  For  the  wave  form  s(t)  which  takes  the  value  of  +1  or  -1  according 
to  s(i)  at  time  (i<=l)d4t  <id,  the  autocorrelation  function  J^(t)  is, 

£(T)-:j=js(t)s(t+T)dt  T=Md 

0 

ForT=0  or  kd     (k;  integer) 

5(kd)=  -^  5"s(t)s(ttkd)dt  =1^  ^dt  I  a, a,,^  =  -  J^ 

If  T=kdtJ^and  0<S'<d,   then  first  (d-^)  of  each  digit  of  s(t) 

corresponds  with  the  last  part  (d-S")  of  s(t  +  kd),  and  the  remaining  ^  of 

s(t)  corresponds  with  the  first  part  (T  of  s(t+(k  +  l)d). 

=  1  I  A 

~   M 
^(T)  is  an  even  function  and  we  can  get  the  same  result  for  any  negative 

value  of Xk  and  the  autocorrelation  function  of  s(t)  has  the  shape  shown 

in  Fig,  2-2. 
ii 


I  ij  I  I  I l_l_J I 1—1 L_! I— L..r  ,  1  J UJ__Lh.C— 1— t 1—1 l—l—! 1 L_l 1 l__i_J I L_J L_ t_J ^— t 


Fig.  2-2  Autocorrelation  Function  of  m-sequence 
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The  power  spectral  density  of  s(t),  G^(^)  is  derived  l?y  using  the 


Wiener-Khintchine  relation. 


0^^=1^8(0  e       dt 


nfO 


I    I    I    I 


ll 


1-1-.. 


-2V 


^ 


£t> 


Fig,  2-3  Power  Spectral  Density  G^(w)  of  m-sequence 
When  we  consider  the  necessary  bandwidth  of  the  band  pass  filter  for 
the  modulated  input  (/2Acos(a;ct+9s(t)),  where  0<6<^  ,  we  can  analjrze  s(t) 
into  frequency  components.  If  we  pick  up  one  frequency,  then  that 
frequency  component  will  be 

M 

Since  M  is  a  large  number,  i^^jk<^^  and  the  amplitude  of  each  frequency 
component  is  supposed  to  be  very  small.  A  phase  modulated  carrier  cOc  by 
a  modu3^ting  frequency  tU^,  with  small  modulation  index  /!?<^1,  is  expressed 
as 

cos(t<;ctt^sinwMt)«  cos  M(t-  /Ssinu^mtsinM^c-t 
This  has  frequency  components  of  cc)^,,  cO(±iU^»     The  frequency  components  of 
the  signal  \/2Acos(u)c.t+0s(t))  are  derived  in  the  same  way  as  for  small 
modulation  index  /S^J,  therefore  in  order  to  preserve  almost  all  power 
of  this  incoming  signal  we  can  pass  only  the  frequency  between  «^c.^T^t 
because  almost  all  the  power  of  s(t)  is  included  inside  of  i^^* 
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CHAPTER  III 
PHASE  LOCKED  LOOP  FOR  CARRIER 

3.1  Coherent  Demodulator 

Since  s(t)=±l,  i^e  can  separate  the  input  signal  ^fZAcosiu^ct^Osit)) 
into  two  parts. 

V^Acos(M;ct+0s(t))=/2Acos6  cosu)^t-i-ii{i)\l?.hsinBr.XnWJi 
If  we  can  lock  the  phase  of  the  local  oscillator  to  the, carrier  frequency, 
we  can  get  the  modulating  signal  '/2Asin0s(t).  While  \fZAcos9cos  cO^t   can  be 
used  for  that  phase  lock.  Hence,  there  will  be  an  optimum  value  of  9  to 
get  the  maximum  signal  to  noise  ratio  in  the  output. 

Phase  jitter  in  the  local  oscillator,  caused  by  the  output  of 
\/2As(t)sinO,  plus  input  noise  causes  the  amplitude  variation  of  signal 
output  in  the  form  of  >/2As(t)sin(9+  9o)   and  acts  to  increase  the  noise 
input  for  the  next  stage. 

The  cross  correlation  function  of  cos(uJtt+^s(t))  and  cos«4t  is 
f(r)=J=l  cos(M^ct+9s(t))cosa;£(t-hT)dt 

~2t1  i[<^os0cos(2uJtt  +  w<:r)-s(t)sinOsin(2a;ct  +  ^cT)  cosuJcXoos^ 

+  s  ( t )  sin6sin  ft^cT  J  dt 
=    fcos^coso^tT-prsin^sintt'cT 
=    icos6cos(a;cT+^) 
<^=  t3.n'(^tsnO) 
<p   is  negligible  when  M  is  large  and  tan^  is  not  excessive  and  we  can 
make  a  phase  lock  loop  like  Fig.  3-1 • 

As  far  as  this  phase-lock  loop  is  concerned,  the  first  low  pass 
filter  can  be  neglected  and  letting  the  low  pass  filter  be  denoted  as 
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V2Acos(aJctf^s(t))/o: 

+n(t; 


l?M 

B.W.  ^% 

?^ 

L.P.F. 

A(JKMS^x;sin(7 

V 

i2Csin 

u)it 

v.c.o. 

li*l   .1'  .    TiV*'/ 

Fig*  3-*l    Phase  Lock  Loop  for  Coherent  Demodulation 
F,("^),  we  can  simplify  this  loop  like  Fig.  3-2. 

V2Acos(a>cttXt  +  0s(t))+n(t)  ACk^cos^sin(Xr'^)+ACk^s(t)sin0cos(Xi- Xo) 

® 


-^Csin(^ct-''X*) 


-n(t)k/2Csin(wct  -t  Ko) 


V.C.O. 


F,{u)) 


Fig.  3-2  Phase  Lock  Loop 
Noise  input  is  assumed  to  be  band  limited  white  Gaussian  noise,  and 
using  narrow  band  approximation 

n  (t )  =  n  (,  (t )  cos  tt>tt +•  n  ^(  t )  sin  uJct 
where  n^(t)  and  na(t)  are  normally  distributed  time  functions  varying 
slowly  compared  withoic*  and  av.(ni)=  av.  (nc)=0,  the  variances 
6'^=^^=H  =  '^\I  W  bandwidth  (cycles/sec) 

•7  one  sided  power  spectral  density  of 
white  Gaussian  noise  (v:attB  sec/cycle) 
When  we  look  at  the  vector  of  this  input  at  the  vicinity  of  the  phase 
locked  point  in  Fig.  3-3 1  \/2Asin0  and  n*  are  found  as  the  causes  of 
phase  deviation,  and  since  they  are  independent,  their  effects  will  be 
evaluated  in  the  sum  of  variances. 


Fig.  3-.;'  Vector  Diagram  for  Phase  Noise 
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By   approximating 
sin(Ai~  Ao)=  \[~  ^o 

in  the  vicinity  of  the  phase  locked  point,  we  can  linearize  the  loop  as 
shown  in  Fig.  3-^«  The  output  of  the  voltage  controlled  oscillator  is 
the  derivative  of  Ao  and  proportional  to  the  input  level.  Thus  we  can 
also  modify  them  in  Fig,  3-^» 


^ 


U) 


AiM4-e5M+^  ^ 


Fig.  3-^     Linearized  Phase  Lock  Loop 


Ao(<^) 


Fig.  3-5  Linearized  Phase  Lock  Loop 
Let       k^Cg.=  K 
Acos0=A, 
then  we  can  get  the  loop  of  Fig.  3-5 •  If  we  neglect  the  effect  of 
0S(tO),  the  solution  for  optirauiii  performance  is  shown  in  the  paper  of 
R.  Jaffe  and  E.  Rechtin.  [5] 

Using  that  result,  F,(aJ)  is  defined  as  F ,iu))=  ■^^—^■^^   and 
the  closed  transfer  function  H(a))  is 


(M 


The  network  for  F,  («^)  is  approximately  constructed  by  RC  components. 
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— o 


RrR.i^->)»^' 


Fron  the  linearized  transfer  function  H(w) ,  we  can  also  derive  the 
error    6M=A{^'v)— /\o^"^) 

=  /^^^-^  B-.^B,%^u.e  (3-1) 

For  the  input  of  frequency  deviation  ^to 

^^^"^^    B/ -e^"k/w  -  u;*  (3-2) 

€(t)=  ^_S^£(.)^'^'^^u)=^^^^^f  t  t>  0 

=  ^  t<0  (3-3) 

Max    £  (t)  is  derived  by  putting  ^=0. 

«(t)L«=  sSfg.  0-6) 


The  effective  bandwidth  of  E(to)   is 


=  ^B,=  1.06Bo    (cycles/sec)  (3-7) 

The  phase  error  due  to  input  noise  is  normal  with  zero  mean  and  the 
variance  is 

3.2  Noise  In  the  Demodulator  Output 

As  shown  in  the  previous  section,  the  output  near  the  phase  locked 
point  is 

Acos9sin6|,+As(t)siJi0cos0„+n^(t)  (3-9) 
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var(n^(t))  =  7^ 

If  we  subtract  the  required  signal  As(t)sin9,  then  the  remaining 
is  unnecessary  and  is  considered  as  noise,  which  is 

Acos0sin^c)+  As  ( t )  sin  6(l"cos^o)^no  (t ) 
a)  The  distribution  of  do 

The  phase  error  Q,  in  eq,(3-9)  is  the  sum  of  the  contribution  from 
the  noise  input  and  the  self  noise  of  s(t). 

Since  H(w)  is  assumed  to  work  in  the  linear  region  and  n^.  is 
normal,  the  phase  error  produced  by  the  noise  input  0„  is  approximately 
distributed  normally. 

The  distribution  of  phase  error  by  the  self  noise  9%  is  derived 
as  follows. 

Consider  a  randomness  property  of  an  m-sequence,  "R-2"  (the  Run 
Property),  "Among  the  runs  of  ONE's  and  of  ZERO'S  in  each  period,  the 
number  of  runs  of  length  i  is  l/Z*"  of  the  runs  of  each  kind,  and  n  = 
log^(M+l)  is  the  maximum  length".  VJhen  we  compare  these  distributions, 
we  can  see  each  digit  has  the  probability  of  |-  to  be  +1  and  i  to  be  -1, 
If  we  take  a  large  number  of  digits  m,  then  the  probability  it  contains 
k  of  +1  is  binomial  distribution. 

P(k  of  +1  out  of  m  digits)  =  (^)(x)^ 
Therefore  if  we  take  m  digits  and  count  the  difference  of  the  number  of 
-1  from  the  number  of  -fl,  the  probability  of  the  difference  to  be  i  is 

P(dif  f  erence  -  i)=  i^  jlzT 

(m  and  i  should  be  both  even,  or  odd,  but  since  m  is  large,  for  a  small 

number  of  i,  we  can  adjust  ra  as  m  or  m-t-1  according  to  the  number  of  i 

for  approximation).  For  a  large  number  of  m,  we  can  approximate  it 

with  normal  distribution. 

Since  a  low  pass  filter  is,  in  some  sense,  an  averager,  we  can 
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conclude  that  ^  is  also  approximately  normally  distributed  with  zero 
mean  and  with  variance  equal  to0Bod»(See  Appendix  l)/^ 

6,,  and  ^  have  different  kinds  of  input  and  they  are  independent, 
therefore,  the  probability  distribution  function  of  ^o'&^^^^  is  also 
normal  with  zero  mean  and  variance  ^^=Se^^+6"g  . 

b)  The  probability  distribution  of  sin^©,  and  (l-cos  ft,)  where  ^^  has  a 
normal  distribution  function. 

Considering  the  noise  output  by  qj,  if  we  want  to  get  good  per- 
formance in  the  next  stage ,  do  should  be  desired  to  be  a  very  small  value 
in  the  variance,  if  not,  this  system  would  not  work  satisfactorily. 
If  do  has  a  small  variance ,  the  contribution  of  60  to  Acos^sin6i)  is 
approximately  evaluated  mainly  by  6b  neglecting  the  other  terms  of  series 
expansion  of  sin^^,  and  for  (l-cos  9o)   by  ^c/z, 

Acos0sin^=  A^cos0 

As  (t  )sin0(l-cos  0o)=  f*As  (t)sin^ 
The  first  one  has  normal  distribution  and  the  second  one  -^  is  Chi- 
square  distribution  with  decree  of  freedom  1. 

Looking  at  Oo   which  is  the  output  of  H(a^),  we  can  see  that  it 

varies  very  slowly  compared  with  each  digit  interval  d,  but  since  the 

ratio  of  maximum  frequency  deviation  of  the  input  to  ^c  is  the  same  as  the 

ratio  of  the  deviation  of  digit  duration  to  d,  and  the  filter  parameter 

Bfl  is  wider  than  the  filter  parameter  tj^  in  the  next  stage ,  Oo  varies  more 

rapidly  than  the  fluctuation  of  phase  in  the  next  stage.  Therefore  the 

effect  of  Oo  to  the  next  stage  is  considered  as  if  Acos^sin6t  were  added 

to  n»  as  random  Gaussian  noise  in  the  approximation  of  (Acos0)9,,  and  since 

AsinO(l-cos  9o)   is  approximated  by  %,,   at  the  output  of  the  next  stage  th^ 

have  the  same  effect  as  the  decrease  of  amplitude  by  a  factor  like  an 

average  of  many  digits, 
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The  averaged  power  of  As(t)sinucos^<j  is 


/ 


:^< 


= ^f^  cr^-"''^'*-  i:w^-^»-'''"'  ^^.j 


=A^sin^0[i^K''*'/^J  (3-13) 

or      (l-cos0^)Asin0^-^^%  Asin0 

[^J  is  distributed  in  %,    distribution,  and  if  we  take  sum  of  n 
digits,  it  becomes  /v„ 

E(^')=l.  var.(^')=2/n 
And  as  n  becomes  large,  the  distribution  of  the  average  becomes  very- 
close  to  1 , 

av.  (/  in0  cos0o)  =  (1-  ^'')Asln9  (3-1^) 

Comparing  this  with  J  7+  Ve  ^  ,  f  oif  sraall  values  of  6^^ 

The  result  shows  good  agreement  \'itii  l;'  .  (3-1^). 

c)  Correlation  of  0©  and  no(t) 

From  the  argument  so  far,  we  could  derive  the  approximate  output  of 

the  di;?cri  in.' t'lr ,  and  It  is  shoAU  as  follows. 

J  td+e""^^^)  (Asina}3(t;  as  the  signal 

-f-(Acosd)9o+  nJt)  as  noise       (3-15) 

What  is  the  distribution  of  noise  as  a  whole? 

Qj  is  uncorrelated  v;ith  cl,  and  so  with  nj>(t). 

Therefore  the  correlation  of  l7(,  and  n«(t)  is  the  correlation  of  q,  and 

nj[t). 

In   a  linear  system  of  transfer  function  E(i*J) , 

h(t)  =  N/2B^e'^*^  cos^  t         for  t>0  (  h(t>=0  for  t<0) 
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6l,(t)=  [^  n(t-T)h('r)dr 
Correlation  of  9»,  and  n|_(t)  is 

r^„a)  =  liB^f  n(t+A)  e„(t)dt 

=  Idjn-ir  (  n(tA)  \  n(t-T)h(r)drdt 

T*»  ^  <    ~T  ■'a 

=  lijti:^  Ph(r)  rn(tU)n(t-r)dtdr 
=  Ch(r)lira ^  \\(t  ^-\ )n(t- T )dtdr 
=  i^h(T)R„a-T)dT  =  j1i(T)7  (^\-t)dZ: 

=A%^^'*'     -4^  (3-16) 

for  X=  0 

>      ■^*'.('')=&  (3-17) 

The  variances  of  ^©and  n«{t)  are 

6;^'=  r^  (3-18) 

The  normalized  correlation  coefficient  f is 


y 


^*'''  S  (3-19) 


Letting  O^AcosS  =  x 

ng.(t)=    y 
E(x)=  E(y)=0 
Ci^'sA^cos^e    6^^ 

6-;=  G-,; 

Changing  variables  x+y  =  u 

y=v 


where  z  =  OTTS^  +■=??- nr+-  ii-        (>c>v      +   ^y^  > 

P(U)  =-==-p=4===:=     .^.ajp    r_  U.^ 1 
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u  =  no(t) 
The  output  is  now  simpliTied  as 


JHU^T^')  (Asin^^(t)-t-  no(t) 
np(t)  is  normal  with  zero  mean  and  variance  5"„, 

The  power  spectral  density  of  this  noise  is  not  flat  in  the  bandwidth 
of  l/d,  because  c?©  has  the  approximate  bandwidth  of  B^.  Let  the  power 
spectral  density  of  n^Ct)  be  G„(«J) 

-CD 

If  we  divide  this  by  Gni'^A         »  then  we  can  get  the  effective  band- 
width  B£,.  The  low  pass  filter  used  in  the  next  stage  has  narrower  band- 
width  than  B^,  then  we  can  use  (~^)  as  the  approximation  of  G^Cto)  because 
the  value  of  Q„{i^)   inside  of  that  bandwidth  o^o  mainly  contributes  to  the 
output . 

Let  the  Fourier  transform  of  n^Ct)  be  HAto),   then  the  Fourier  transform 
of  fAcos^la  is   ^*^(^>[acos0)  0;w) 

(acos9)@»(^^=  (Acos^l^iO 

GhM  =  [l%H^■#-'H(u;)|7  =  2.917 
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CHAPTER  IV 
DEUY  LOCK  DISCRIMINATOR  FOR  A  BINARY  SEQUENCE 


4,1  Tracking  and  Acquisition  Performance  of  a  Delay-lock  Discriminator 
This  has  been  solved  by  J.  J.  Spilker,  [8 J  and  here  we  would  like  to 

follow  the  general  results.  This  section  explains  the  tracking  and 

acquisition  performance,  and  the  next  section  shows  the  effect  of  noise 

on  the  performance. 

If  we  take  the  difference  of  the  outputs}  from  two  crosscorrelators 

which  differ  two  digit  intervals  in  Fig.  4-1,  and  separate  it  into  an 

averaged  value  A^kDjCe^)  and  another  fluctuating  part  (self -noise) 

A(jkn^(t,6),  these  are  derived  as  follows. 

Input         ^^^  Loop  Filter 

Aos(t-HT)fno(t) 


ks(t+T  -d) 


gfF^iv) 


^ 


M-l 


ks(t+T+d) 


«-2  -  -  —  — 


Voltage 
Controlled 
Clock 
gain  gc 


n-stage  Feedback 
Shift  Register 

Fig.  4-1  Delay  Lock  Discriminator 

The  autocorrelation  function  of  an  ra-sequence  is  already  shown  in 

Fig.  2-1,  and  from  that  figure  we  can  draw  D^(6)  as  this. 
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D^(6)=  E[s(t+T)s(t4-dfT)-s(t+T)s(t-d-HT) 

=  {R,(f-d)-R^(<f+d)]  (4-1) 

The  self -noise  in  the  output  of  the  multiplier  network  is 
n^(t,  G  )=  s(t-<-d^-$)s(t-l-T)-s(t-d+T)s(tH')-D^(6)     (4-2) 
For  f=kd,  k:  integer,  using  the  delay-and-add  property, 

a)  k=0,  (modulo  M) 

n^(t  .kci)-  s(t-Hdi-T).s(t+T)  -s(t-d-t-T)s(t+T) 

=s(tfT+d)s(ttT)-s(t+-T-d+d)s(ttT-d) 

=s(t-»-T+jd)-s(t+-T-d-hJd) 

=  s(tt'B3d)-s(t-»-Tt(j-l)d)  (4-3) 

b)  k=+l,  (modulo  M) 

n^(t,kd)=s(t+d^-T+d)s(ttT)-s(t-dtTtd)s(t+T)-D^(id) 
=  |l-s(t-hTtid)-4tl_ 

[s(t-HT-Md)-l+-^ 
=i[s(tmid)--^J  (4-4) 

c)  k= other  integer 

n^(t.kd)=s(t^H-rd)-s(t-HrHd)  (4-5) 

Further  discussions  of  this  self-noite,  i.e.  its  power  spectral 
density  and  the  output  from  the  loop  will  be  explained  in  the  next  section. 
A)  The  output  of  the  crosscorrelator  network  is 

kAo[DM(t)fn^(t,6)fn«(t)/Aj 
This  is  filtered  by  a  low  pass  filter  gfFzit*^)   and  drives  the  voltage 
controlled  clock  generator.  This  relation  is  now  shown  in  the  following 
Fig.  4-2.  Obviously  the  control  part  of  the  crosscorrelator  output  D^(€-) 

A 

is  a  function  of  T-T=^,  and  T  is  the  integral  of  the  controlled  frequency. 


T  =  ^kA,gcgfFa('<^)[D^(6)fn^(^,6)fno(w)/Ao] 


28 


AoS(t-»-T)^-no(t) 


■^   Crosscorrelator 


kAoDa(e>Mi^<-no  /a. 


gfF^(«t>) 


ks(t+-T) 


VCO 


Fig.  4-2  Simplified  Delay  Lock  Discriminator 
Li  the  locked-on  state,  €<d, 


M 


€=T-T 


0^=dkA.g,g,^F»ri  +  |^;±^] 


^fl^   -h/TMJoJtO 


Using  the  same  kind  of  filter  as  used  in  Chapter  3. 

go=  kAegtgf 

(ft  -a**  L  £<  (Mtl'>/M         J 

The  closed  loop  function  Eito)   is 


^ 


The  value  of  T/d  without  noise  is 

T 


Letting 


f(t) 


=    at 


w' 
=  at-^e"^^  sin^t 


'vf£ 


-i^t 


iC^^^^^^^J^ir 


w, 


max 


gg). 


at  t  = 


TL 


(4-6) 


(4-7) 

(4-8) 
(4-9) 


(4-10) 


(4-11) 
(4-12) 
(4-13) 


The  critical  value  of  m^o  for  expected  maximum  value  of  a  is  to  be 
decided  frcMn  nonlinear  analysis.     If  we  limit  this  as 

max.  <?/d  ^  1 
then  a/wo^2.19  (4-14) 
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B)     Acquisition  transient 


(^-15) 


If  we  take  D^(^)  onlxi  from  the  previous  equation. 


I    +   fzJio/uJo 


L^-TuT 


Djfe)- 


M 


A71-/ 


(4-16) 


Let 


iS    implies  "i^   =^ 
Da(£)7Crff  is  denoted  D(x) 


g^j  A/*  ==  g 


D(x)=J  2-x 

X 

-2-x 

io 

f  =D.(x)=|    1 

-1 

.   0 


1^  x^  2 
-1  ^  x  ^  1 
-2^x  <-l 
otherwise 
-1^x41 
-2^x^-1 
otherwise 


or 


l^x<2 


then  (y^^|^)(r^)=('  +  /^||f)Dr^) 
y  +  y/g  =  X -I- x/g -t- D(x)  f  V2D' (x)x 


(4-17) 


where  x  =  44 


etc. 


(4-18) 


(4-19) 


x  +  (j+N/2D'(x))x+-D(x)=y^y/g 
The  slope  of  phase  trajectory  T  is 

If  we  plot  this  performance  on  the  phase  plane  x  and  x,  D* (x)  and 

D(x)  are  the  function  of  only  x  and  we  can  find  the  dividing  line  on 

x-=-2,  x=-l,  x  =  l  and  x— 2.  Let  these  regions  be  called  A,  B,  C,  D,  E 
as  shown  in  Fig.  4-3. 
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X 

(K 

B 

c 

c 

D 

& 

~x 

-1 

o 

1 

2. 

Fig.  4-3  Regions  in  Phase  Trajectory 
Region  A  and  £ 

D'(x)=D(x)=0 


x  +  gx=y  +  |y  ^ 

r4  '/^ 


(4-20) 


IT  does  not  contain  x,  but  is  defined  only  by  x  and  the  input. 
Thus,  if  x  =  y  at  the  initial  point,  then 


» -^1^ 


X 

y  =  x 

x-y+-Xo         x^:  initial  position         (4-21) 
If  X:f  y  at  T=0,  letting  the  initial  values  of  x,  x,  y,  y  be  x^,  x»,  0, 
y^  ,  and  substitute  them  into  the  general  solution  of  this  equation 

x  =  c,-i-Cje  <'+  y 
and  we  obtain 

x  =  y->-x,+  g(xo  -y^)(l_e  ^)  (4-22) 

Region  B 

D'  (x)=  -1 

D(x)=  -2-x 

For  a  constant  velocity  input,  i.e.  y  =  0,  the  point  having  a  constant 
T«  k  draws  a  line . 


(ki-^-  ?)x  -(xi-2-^y/g)  =  0 
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If  g  is  large  enough  to  approximate  l/g=0,  then  upper  part  of  the  line 
combining  (x= -1,  x=-0.518),  (x=-2,  x=0)  converges.  This  line  is 
derived  as  follows, 


Fig,  ^-4  Phase  Trajectory  in  Region  B 
for  constant  T=k 

(k-\/2)x  -  (x  +  2)  =  0 


i-r- 


Zi-X 


\fz-k 


I 


X'-j 


^^2-k 


k<0 


the  slope  of  the  line=  k  =  -^  =-  '  . 


therefore 


k^  -{2  k  -1*0 


Region  C 

D(x)=  X 

D'(x)=l 

0  ~  -^ 


(if-24) 


X  = 


^  ^  m  ±ji. 


Region  D 

D(x)=2-x 
D'  (x)=  -1 


(^-25) 


r  -t-  '4  -  v/i: 

The  line  which  divides  the  plane  into  converging  and  diverging  parts  is 
derived  by  the  same  way  as  that  in  i-egion  B, 
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Fig.  ^-5  Phase  Trajectory  in  Region  D 
The  result  obtained  by  CDC  l60^  digital  computer  i>r  ^  =  03   end 

several  initial  conditions  is  shown  in  Fig,  ^-6,  (Coniputer  programming, 

see  Appendix  II). 

From  this  figure  x  =  2,35»  x=0  at  the  initial  condition  is  the  limit 

to  converge.  It  implies  that  -^<2.35» 

Time  response  of  position  error  and  velocity  error  are  shown  in 

Fig.  4-7,  8.  From  these  figures  the  locking  time  is  approximately 

The  phase  trajectory,  x  vs. t and  x  vs.Tfor  g-lO  are  shown  in  Fig, 
4-9 t  10,  11.  From  these  figures  we  can  find  the  position  error  of  0.2, 
and  the  limit  of  y  is  approximately  2.0,  Lock  time  is  also  S/u^ti^Zir /uJo  . 
4,2  Noise  and  its  Effect  on  the  tracking 

In  this  section,  I  would  first  like  to  discuss  the  effect  of  self- 
noise;  second  the  effect  of  the  input  Gaussian  noise,  and  finally,  the 
effect  of  the  limiter  before  this  stage, 
a)  Self-noise 

Self -noise  is  a  function  of  time  and  phase  error.  Continuing  from 
the  last  equation  (4-2),  power  spectral  density  for  the  time  error  ^ 
being  an  integer  multiple  of  d  is  derived  as  foiiovjs; 

n^(t,kd),  where  k=0 

The  autocorrelation  function  of  this  self -noise  is  the  sum  of  three 
triangular  waves  as  shovm.  Ely  the  Wiener-Khinchin  relation  the  power 


33 


-^ 


spectral  density  is  the  Fourieij 
transform  of  this  wave» 


G„^(w,f=0) 


=  2G^(u>)(l-cos««)d) 


(^-26) 


Since  n^(tt,d)  is  the  same  as  signal « 

ksr  other  integer 

By   transforming  the  autocorrelation  function 

G„^  (u)„6=kd)=  2G^(«;)[l^cos(r=q)u>d]  (4=28) 

We  are  interested  in  the  effect  of  this  noise  in  the  linearized  region 

-d^e^d. 

For  the  extreme  case  6=±d^   the  power  of  self=noise  in  the  output  of 

H(u>)  is 

-d^OJod  (it-29) 

For  the  least  extreme  case„  (£  =o 


b)  The  mean=>square  delay  error  caused  by  the  input  noise 

4,3  The  effect  of  input  limiting 

By   the  previous  discussion  we  approximated  the  input  noise  as 

2 

Gaussian  with  zero  mean  and  the  variance  6'„o  The  signal  is  a  square 
wave  of  amplitude  ±Ao» 


Given  a  digit  which  has  equal  probability  of  +A  or  -A  with  additive 
Gaussian  noise,  and  if  we  decide  the  signal  is  ±A  according  to  the 
polarity  of  input,  the  error  probability  for  each  digit  is  derived  as 
follows : 

When  we  defined  A=Ji(l+e  **")  Asin^,  this  was  derived  as  an 
averaged  value  for  some  interval,  hence  if  we  want  to  get  the  probability 
of  error  for  each  digit,  we  can  not  use  A©  as  the  aiiip3j.tude  of  the  signal. 

The  probability  of  error  Pg.  is 

Pe=  P(s(t)=+l)P(output  is  negative/s(t)  =  +l) 

P(s(t)=  -l)P(output  is  positive/s(t)=  -1)     (4-31) 
The  conditional  probability  of  each  term  is  the  same  value,  because 
additive  noise  has  zero  mean  and  the  signal  amplitude  is  the  same  in 
absolute  value. 

PE=2P(s(t)  =  -»-l)P(output  is  negative/s(t)=-»-l)     (4-32) 
The  output  of  the  first  discriminator  is 

Acos^sin^o ■•"  As(t)sin0cos6*+-no(t) 
Qo  and  Kjfl  are  Gaussian,  with  zero  mean  and 

they  are  not  mutually  independent,  but  the  normalized  correlation 
co-efficient  f is 

From  eq.  (4-32) 

P(output  is  negative/s(t)  =+l) 
=P  [Asin(  ^  f  ej  +  nc  <  0  J 
=  P[nc^  -Asin(e-t-a>)] 


2TrS-e(5'^ 


TirW  f-^'-'-  %  )^f  77^(- 1^  '*'^«'-  ^  V^^  ('^-33> 
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where  erf(x)=^)^-^  Ax. 

Consider  a  pulse  train  of  ra  digits,  where  each  digit  has  probability 
of  error  Pg,  then  if  all  pulses  are  1,  the  averoge  of  this  pulse  train  is 
derived  as  follows. 

The  number  of  error  k  in  m  digits  distributes  binormally,  and  if  m 
is  large  enough  and  mPg  is  not  small,  then  this  distribution  can  be 
approximated  by  a  normal  distribution  with  mean=  mP^ ,  variance 
mPpCl-Pg.),  In  this  case  each  error  causes  the  decrease  of  2/ra  in  the 
average,  and  this  decrease  is  normally  distributed. 
The  decrease  of  the  averaged  voltage  v  is,  thus 

If  the  pulse  train  is  composed  with  the  same  number  of  -f-l  and  -1, 
i.e.  )']/?.   oi  positive  pulses  and  ra/2  of  negative  pulses,  then  an  error 
which  occurred  in  positive  pulses  causes  negative  error  2/m  in  the 
average  voltage  and  the  same  for  an  error  in  negative  pulses. 

XI     positive  error  voltage  caused  by  errors  in  negative  pulses 
y:  negative  error  voltage  caused  by  errors  in  positive  pulses 
Since  the  error  for  each  digit  is  independent,  x  and  y  are  independent 
because  they  are  the  error  caused  by  only  negative  or  positive  ^lulsey. 

Since  they  are  independent  „  the  probability  distribution  of  u  =  (x+y)  is 
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This  variance  is  equal  to  that  in  the  positive  pulse  train. 

The  effects  these  digit  errors  cause  in  the  next  binary  discriminator, 
will  be  discussed  in  two  phases,  first  at  the  locked  point,  and  this  can 
be  considered  as  a  phase  jitter  in  the  phase-locked  position,  and  second 
at  the  limit  where  the  phase  error  becomes  ±d.  This  may  be  considered  as 
the  modification  of  D^(6)  curve.  These  two  effects  will  show  the  analysis 
of  the  system, 

^Vhen  one  digit  error  occurs  in  the  input,  how  does  this  appear  in 
the  outputs  of  two  crosscorrelators?  If  the  system  is  locked  in  phase, 
^  —  0 ,  the  outputs  are  also  a  cyclic  permutation  of  that  m-sequence  and 
their  phases  differ  by  1  digit.  Since  the  autocorrelation  function  of 
this  sequence  is  nearly  zero  for  e=i:d,  the  errors  introduced  in  the 
outputs  of  the  crosscorrelators  do  not  have  any  relation  to  their  polarity 
except  the  fact  that  they  occurred  at  the  same  time.  About  the  polarity, 
from  the  characteristics  of  the  randomness  of  m-sequence,  the  probability 
that  both  errors  are  positive  or  negative  is  ^  each,  and  the  probability 
the  sum  is  zero  is  j.     This  implies  the  error  probability  is  decreased 
by  a  factor  of  2,  and  we  can  use  ?-e/2   instead  of  Pg.,  and  the  change  in 
the  average  voltage  output  caused  by  one  error  becomes  double. 

Thus,  we  can  modify  the  equation  of  the  noise  voltage  as  follows. 

For  the  case  of  6^d,  one  error  always  appears  as  negative,  and 
another  as  positive  of  negative.  This  time  we  do  not  need  to  modify 
Pg ,  but  we  can  add  two  effects. 


We  have  been  discussing  the  effect  of  noise  on  the  average  through 


m  digits  so  far,  and  now  we  have  to  make  sure  of  the  meaning  of  the 

average, 

^3    'K 


The  inverse  Fourier  trfnsforin  of  H(^),  h(t)  is 


-7^ 


tr 


h(t)  =  N[2a;„e  "^  cos^'t 
The  response  of  this  filter  to  a  pulse  of  height  1  and  duration  d  is, 
if  we  take  the  response  after  time  d, 

^e'%sin(^t^f-) 
Since  the  chance  of  error  is  very  small,  it  will  be  improbable  that  the 
digit  in  error  will  be  discovered.  Therefore  we  can  neglect  the  response 
of  this  filter  in  the  digit  duration  and  assume  this  as  an  impulse  response. 

When  we  look  backward  for  a  time  interval  of  T,  standing  at  an  arbi- 
trary time  t|,  if  we  use  a  real  averaging  device,  the  value  at  tjis, 
assuming  there  are  k  impulses  to  occur  in  T, 

u,=Y[Zd  J^Ct-t:  )dt=4i^  . 
If  we  use  H(^)  as  an  averager,  the  output  average  differs  depending  upon 
the  positions  of  impulses,  but  when  we  think  of  the  interval  of  each 
impulse,  this  is  approximately  distributed  exponentially,  and  after  taking 
many  intervals  T  containing  k  impulses,  we  find  k  impulses  are  distributed 
in  equal  intervals  as  an  average. 

Thus,  the  output  average  is  the  integration  of  that  response.  The  time 
internal  T  should  be  determined  for  H(^),  and  looking  at  the  response 


'^e  ^  N/2'sin(gtt^),  we  can  neglect  the  effect  of  the  tail  after  t=f^. 


T-^f  (^-37) 

u.^^^^^'e^  Nf2sin(gt^^)dt 

-   3w^     ~^4    T  (^-38) 


Thus  we  can  approximate  the  effect  of  H(tt))  by  the  effect  of  averaging. 


CHAPTER  V 
THRESHOLD  FOR  THIS  SYSTEM 

5,1  Summary  of  the  System  Performance 

When  we  evaluate  the  performance  of  a  phase  lock  system,  maximum 
frequency  deviation,  or  in  other  words,  maximum  target  velocity  that  a 
tracking  system  can  acquire,  and  transient  time  for  phase  lock,  probabili- 
ty of  cycles  slipping  and  r.m.s.  phase  error  at  the  steady  state  will  be 
considered.  The  problem  of  the  cycles  slipping  is  considered  very  im- 
portant for  this  system,  because  the  acquisition  time  in  A,  E  region  in 
Fig,  4-3  can  be  very  large.  We  can  not  allow  cycle  slip  in  the  delay 
lock  discriminator.  Cycle  slip  is  related  to  r.m.s,  phase  error  at  the 
steady  state.  From  the  experimental  results  of  V/,  F.  Roland,  [7l  we  can 
find  the  linearity  of  6"^  vs,  B©  up  to  20-30  degree,  and  the  expected 


2.TT 

number  of  cycle  slip  per  -^  is 

6.10"^  for  6*9.  :  20" 

6.10'"^  for  6-©, ;  25'' 

2,35il0"^         for    6"^^  ;  30° 

— ^t  r~ 

If  we  use  equation  (3-3),  the  time  for  e  ^  =0.1  is  ^'^i    '^^   . 

This  is  for  the  initial  position  error  of  zero,  and  we  can  approximately 

use  twice  the  value  of  it,  ^  ,  as  an  acquisition  time  frcM  the  initial 

error  of  ic  , 

Let  the  expected  number  of  cycles  slip  per  ^^/r   be  p^  and  the 

settling  time  be  "^^/Bq  ,  then  as  an  average 

is  the  number  of  cycles  slip  in  one  period.  If  one  cycle  of  slip  occurs, 
-g-  in  M  is  in  transient  responce  and  it  can  cause  longer  transient 

^5 


period  of  the  next  stage. 

For  the  cycles  slip  in  the  next  stage,  Spilker's  paper  shows  that 
for  delay  error  having  6),^0,30d,  a  probability  to  lose  the  locked-on 
state  in  the  absence  of  transient  errors  is  negligible  from  experiments. 

While  cycles  slip  determines  the  limitation  for  the  values  of  the 
s/N  ratio ,  d  and  the  bandwidth  parameter  Bo ,  k)o  ,   the  maximum  frequency 
deviation  restricts  the  lower  limit  of  B©  and  ^4.  From  the  phase 
trajectory  of  the  transient  response  of  a  phase  lock  loop  for  a  sinusoid 
by  A.  Viterbi ,  the  maximum  frequency  deviation  at  position  error  =  0  to 
be  acquired  is  approximately  2.I5B0  for  this  phase  lock  system  concerned. 
If  we  restrict  ^(fc)l<l  in  the  linearized  system. 


'rtlfly 


^^2.19  from  eq.    (3-6) 

or  ^<2.15  (5-2) 


B, 


o 


For  the  delay  lock  discriminator,  from  Fig.  k-6   and  4-10, 
i^  2 
This  is  §^^  2.  (5-3) 

Or,  using  the  result  of  eq.  (4-14), 

i^  2.19  (5-4) 

The  transient  time  for  phase  lock  is  decided  mainly  1:^  B^  and  <t)^ , 

as  follows, 

2.Tr 
"^lock  ~  ~W  ^^^  ^^^   phase  lock  loop 

These  transient  interval  should  be  limited  by  the  actual  situation. 
5.2  An  Example  of  the  Calculation  Procedure 

We  are  going  to  decide  0,  Bo  and  cOo  for  a  radar  using  this  phase 
lock  system.  The  specification  of  this  radar  is: 

Range  500  Miles 


46 


Carrier  frequency         200  Megacycles 

Pulse  duration,  d         5  Microsecond   (25OO  feet) 

Number  of  digits  in  one  period    M  =^1200 

T=Md  =  5  Milisecond 
Maximum  target  velocity,  ir  3  I'lach,  '}0Q0   feet/sec 
Liput  s/N  ratio  less  than  0  dB 

a)  Lower  limit  of  B^, 

Bo^^  =^00  rad/sec 
for  this  Bo,  tj^  .=  -^  =  1.6  msec 

This  will  be  short  enough  couiparing  with  a  period  of  Md-'5  msec  for 
the  first  acquisition  sjid  we  do  not  need  to  increase  B,>« 

For  the  delay  lock  discriminator  we  put 

T(t)  =  adt  in  eq.  (4-11), 

_  2V  I  J.  /  5 

4^2.0 

"^o^  0.6  rad/sec 
but  t,,,^=^  i5'5) 

for  fJo-0,6        t,  .  =  10  sec 

This  is  too  large  for  this  operation,  and  we  need  to  make  tJo  large 
in  the  limit  where  the  cycle  slip  does  not  occur. 

b)  The  number  of  cycle  slips  per  ^ 

B,=  4000       d  =5.10       M-1200 

Once  a  cycle  slip  occurs,  it  takes  1.6  msec  to  recover  and  this 
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becomes  a  loss  of  about  one  third  of  a  period.  Thus  we  want  to  limit 

1^  in  3.10  ,  or  r.m.s.  phase  error  up  to  25*. 

c)  ^%=^e„-^^ff^ 

f-^^J. ,'^/n  --       '         If  ''"^^^      ) 

where  K  is  input  s/n  ratio 

6->B,d0' 
6o  61  and  ^  vs.  Q    for  various  values  of  K  are  shown  in  Fig.  5-1. 
d)-^^i(l^e~^^'^')sin^d 

Under  the  limitation  of  6g^<0.19  (  =  (25°f  ),  we  can  plot  4l  vs.  9 
for  various  K,  as  shown  in  Fig.  5-2. 

From  Fig.  5-2 1  we  can  find  9  which  maximizes  ^^/^^  with  respect  to 
K,  and  this  is  in  Fig.  5-3 •  The  maximized  Ae/A*"  vs.  K  is  shown  in  Fig.  5-^« 
e)  The  mean  sequare  delay  error  ^„  is  the  sura  of  the  contribution  of 
noise  and  self -noise.  Taking  the  worst  value  of  self -noise,  from  eq.  C^- 
30)  and  (4-29), 

In  Fig.  5-5 1  we  showed   *^/4^  vs.  K  v/ith  inaxirtiurti  value  of  ^ol^  t   for 
various  values  of  ^o.  In  order  to  avoid  cycle  slip,  we  limited  ^/d^O.2, 
and  the  maximum  value  of  ^c»vs  K  to  satisfy  this  limitation,  and  t|„^^  vs. 
K  for  this  tJoare  finally  shown  in  Fig.  5-6. 

From  these  figures,  we  can  see  that  this  system  can  not  show  good 
performance  for  the  S/N  ratio  less  than  -15  dB, 

5.3  The  effect  of  a  limiter  before  the  delay  lock  discriminator  does  not 
improve  s/N  ratio,  because  for  low  s/N  ratio,  %  approaches  to  -f  and  phase 
lock  control  voltage  D^(€-)  effectively  decreases  to  zero. 
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CHAPTER  VI 
CONCLUSION 

The  result  in  Chapter  5  shows  that  if  we  use  the  demodulator  in  the 
limit  to  avoid  cycle  slip,  the  amplitude  of  the  output  signal  A^  is  mainly- 
decided  by  sin  0  ,  and  r.m.s.  phase  error  <5^^  is  also  determined  mainly  by 
the  noise  input  near  the  upper  limit  of  Z^-^y^^  .  Eq.  i^'j)-'^)   shows  that  the 
variance  of  delay  error  ^^^is  the  sum  of  input  noise  and  self -noise  and 
both  are  proportional  to  ^o,  therefore  for  better  s/N  ratio  the  self -noise 
is  dominant.  We  are  interested  more  in  the  low  s/N  ratio  region  than  in 
the  high  one,  and  we  can  also  approximate  6;  by  the  contribution  of  the 
input  noise. 

When  we  want  to  design  this  phase  lock  system,  we  will  find  several 
opposing  factors. 

1)  The  phase  lock  ability  over  the  range  of  the  frequency  deviation 
limits  the  lower  side  of  Bp,  ^o, 

2)  The  cycle  slip  bounds  the  upper  limit  of  ^  and  Bo,  '<Jo  with  respect  to  K. 

3)  For  a  satisfactory  acquisition  time,  we  have  to  get  large  B^and  oio  • 
k)     The  effect  of  d  is  related  to  the  input  bandwidth  and  '^o*  d  is  also 
related  to  the  range  resolution  of  a  radar.  Small  d  will  offer  a  good 
resolution,  but  on  the  other  hand  it  costs  s/N  ratio  and  limits  the  input 
frequency  deviation. 

We  tried  to  solve  the  problem  for  a  constant  s/n  ratio  in  Chapter  5, 
Actually  input  S/N  ratio  can  not  be  constant  in  a  system  application,  and 
we  can  estimate  only  the  least  desirable  value  or  a  range  of  s/N  ratio. 
When  we  use  this  system  for  ranging,  6"^  does  not  have  much  meaning  in 
itself,  but  is  very  important  as  a  cause  of  cycle  slip.  Looking  through 

5^ 


these  four  factors,  2  and  3  are  Especially  significant  in  opposing  each 
other,  and  we  will  finally  decide  Bo,  and  0)0  by  applying  a  weighting 
factor.  The  actual  designing  process  must  be  decided  according  to  each 
case,  and  we  want  to  point  out  the  fundamental  relations  between  the 
above  four  factors,  for  they  must  be  further  investigated. 
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APPENDH  I 
SELF-NOISE 


The  phase  error  output  by  signal  s(t)  itseiLf  is  9s{'<>)Ei'^) 
considering  the  relation  shown  in  Fig,  3-5 • 

Since  M  is  very  large,  ^  is  very  small  and  we  can  approximate  this 
relation  by  an  integral  form. 

^»'     ■•■  —  "JIT  JLco  ^vt  /    \      Bod'i/2.        I 

T  _  iJoci'j.f' __^2lti^ ., 

1=1,-1, -I3 

Moreover,  let    g(z)=^'i  ^^^/^^^^z^ 

and  consider  contour  integrals  along  the  paths  a,  r,  H  or  Z^, 

Since  B^d  is  positive,  the  contributions  of  the  integral  path  H   to  the 

«    /  'a  jJ^^ 

contour  integrals  Cg(z)dz,  and  i  f g(z)e^  *   dz,  and  Pz.  to  i5g(zi)^   dz 
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are  all  zero.  The  integral  along 
the  path  Tdoes  not  converge  to  zero 
for  each  integral  as  t-tQ,   but  the 
sum 


jBc<f'i 


-/^'^ 


converges  as  follov/s, 


liinfyg(z)dz-ijg(z)e^'^'''*dz-i^g(z)e"^'^'^'^dz| 


'0 


Therefore  1  =  1,  -  I^^-  1^ 


IM'& 


=  2Wj  \  ^_  (residue  for  g(z))  -  i  21  (residue  for  g(z)e"'''^''  ) 

"*  UHP  UUP 

-H  i   "2!    (residue  for  g(z)e~^         ,  including  z  =  0) 

l-Hp 

For  I,,  the  sum  of  residues  is  4?j, 


-^^"^hsin^' 


ftc« 


-h     cos^j 


residue  of  z=0  is    -jj^ 


SA' 
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If  Bod  is  small  enough  to  neglect  second  power, 
I  =  l  +  i^fl-3(l-^^)|#  -(l-^^)j 
=^#M  =  1.06B„d 


2. 
'©5 


Thus  ^'■=9^Boci 
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APPENDIX  II 
COFxPUTER  PR0GRAI4  FOR  EQUATION  (4-18) 

PROGRAM  ACQUIRE 

DIMENSION  X(30),XJ:'0T(30),C(15) 

C(10>=1. 

1  CALL  INTEG1(T,X,XD0T.C) 
IF  (X(l)+2.0)  10,11,12 

12  IF  (X(l)+1.0)  11,13,14 

14  IF  (X(l)-l.O)  13,13,15 

15  IF  (X(l)-2.0)  16,16,10 

10  DX=0.0 
DXDOT=0.0 
GO  TO  2 

11  DX=-2.0-X(l) 
GO  TO  3 

16  DX=2.0-X(1) 

3  dxi")Ot=-;l.o 

GO  TO  2 

13  DX=X(1) 
DXDOT=1.0 

2  ACCE=C(1) 
SPED=C(1)*T+C(2) 

ERR0R=C(3)*SPEDfACCE-(C(3)+1.4l4213*DXD0T)*X(2)-DX 
XD0T(1)=X(2) 

XD0T(2)=ERR0R 

GO  TO  1 

END 
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END 


For  g= infinity 

time;     from  0.  to  9*0*  with  interval  0,005 
C(1)=0. 

C(3)=l/g=0. 

Initial  value  of  X(l)--2.0 

X(2)=^y 
For  g=10 

time;  from  0,  to  14.0,  with  interval  0,01 

C(1)=K) 

C(2)=y 

C(3)=l/g=0.1 

Initial  value  of  X(l)=-2.0 

X(2)=y 
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APPENDIX  III 
SUBROUTINE  INTEG  1 

This  subroutine  provides  a  simple  but  moderately  flexible  means 
for  the  solution  of  ordinary  differential  equations,  employing  integra- 
tion by  the  fourth-order  Runge-Kutta  numerical  method.  Input,  and  both 
print  and  graph  output  routines  are  built  in.  The  step  size  of 
integration  can  be  changed  by  appropriate  data  input  and/or  programming, 

Title:  Runge-Kutta  Solution  of  Ordinary  Differential  Equations 
with  Built-in  Input  and  Output  Routines 

CO-OP  ID:  D2-NPS-INTEG1        (FORTRAN  60) 

Programmer:  J.  R.  Ward 

Date:  June  1964 
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